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Abstract. Assuming certain forms of the stream function inverse solutions of an 
incompressible viscoelastic fluid for a porous medium channel in the presence of Hall 
currents are obtained. Expressions for streamlines, velocity components and pressure 
fields are described in each case and are compared with the known viscous and second- 
grade cases. 
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1. Introduction 

In recent years the theoretical study of MHD channel flows has been a subject of great 
interest due to its widespread applications in designing cooling systems with liquid met- 
als, petroleum industry, purification of crude oil, polymer technology, centrifugal sepa- 
ration of matter from fluid, MHD generators, pumps, accelerators and flow meters. The 
results of these investigations are not applicable to the flow of ionized gases. In an ion- 
ized gas where the density is low and/or the magnetic field is very strong, the conductivity 
normal to the free spiraling of electrons and ions about the magnetic lines of force before 
suffering collisions; also, a current is induced in a direction normal to both the electric and 
magnetic fields. The phenomena, well-known in the literature, is called the Hall effect. 
The study of magnetohydrodynamic flows with Hall currents has important engineering 
applications in problems of magnetohydrodynamics generators and of Hall accelerators 
as well as in flight magnetohydrodynamics 1911011 1I18I20I21I23I24I . 

An understanding of the dynamics of fluids in porous media has practical interest in 
such disparate fields as petroleum engineering and ground water hydrology, with applica- 
tions ranging from hydrocarbon migration in reservoirs via packed-bed chemical reactors 
to agricultural drainage and irrigation. In the widely used continuum approach to trans- 
port processes in porous media, the differential equation governing the macroscopic fluid 
motion is based on the experimentally established Darcy's law |2|, which accounts for 
the drag exerted by the porous medium. Brinkman |5| studied Darcy's law by adding to 
it a viscous term in order to account for the vorticity diffusion caused by the boundary 
resistance, whereas the combined influence of inertia and viscous effects on the flow and 
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heat transfer in the vicinity of an impermeable plane surface was discussed by Vafai and 
Tien |28| and Kaviany 1 13 1. Kaviany's analysis was for steady-state flow, the oscillatory 
motion was studied by Khodadadi 1 1 41 131 and transient fluid motion in a porous medium 
channel was discussed by Anderson and Holmedal 1 1 1. 

The governing equations that describe the flow of a Newtonian fluid is the Navier- 
Stokes equations. These equations are non-linear partial differential equations and known 
exact solutions are few in number. Exact solutions are very important not only because 
they are solutions of some fundamental flows but also because they serve as accuracy 
checks for experimental, numerical and asymptotic methods. Since the equations of 
motion of non-Newtonian fluids are more complicated and non-linear than the Navier- 
Stokes equations, so the inverse methods described by Nemenyi 1 17 1 have become attrac- 
tive. In these methods, solutions are found by assuming certain physical or geometrical 
properties of the flow field. Kaloni and Huschilt 1 12 1, Siddiqui and Kaloni |25|, Siddiqui 
1 26 1, Benharbit and Siddiqui 0, Labropulu 1 16 1 and Siddiqui et al\21\ used this method 
to study the flow problems of a second-grade fluid. 

In this work, we discuss the effects of Hall currents on the steady flow of an electrically 
conducting second-grade fluid in a porous medium channel. For such a fluid equations are 
modeled for a grade of fluid two and are solved by assuming certain form of the stream 
function. The graphs are plotted explicitly in the functional form to see the behaviour of 
the flow field. 

The paper is arranged in the following fashion: In §2, governing equations and formu- 
lation of the problem are given. Section 3 consists of two parts. First part is the general- 
ization of Siddiqui's 1 26 1 work and the second part deals with some special flows called as 
Riabounchinsky type flows and finally, concluding remarks are given in §4. Stream func- 
tion, velocity components and the pressure fields are derived in each case. Moreover, the 
streamlines are plotted in each case to see the flow behaviour. 

2. Governing equations 

The constitutive equation of an incompressible fluid of second-grade is of the form 1 22 1 

T = -pI+y,A 1 + a 1 A 2 + a 1 A 2 l , (2.1) 

where T is the Cauchy stress tensor, —pi denotes the indeterminate spherical stress and 
jU,0Si and 0C2 are measurable material constants. They denote, respectively, the viscosity, 
elasticity and cross-viscosity. These material constants can be determined from viscomet- 
ric flows for any real fluid. Ai and A2 are Rivlin-Ericksen tensors |22| and they denote, 
respectively, the rate of strain and acceleration. Ai and A2 are defined by 

A 1 = (gradV) + (gradV) T , (2.2) 
A 2 = ^+Ai(gradV) + (gradV) T A 1 . (2.3) 

Here V is the velocity, grad the gradient operator, T the transpose, and d/df the material 
time derivative. 

The basic equations governing the motion of an incompressible fluid are 

V-V = 0, (2.4) 

p^=JxB + divT-^V, (2.5) 
at K 
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V ■ B = 0, V x B = n m J, V x E = 0, 



(2.6) 



where p is the density, J the current density, B the total magnetic field, fi m the mag- 
netic permeability, E the total electric field current and K the permeability of the porous 
medium. Making reference to Cowling [6 1, when the strength of the magnetic field is very 
large, the generalized Ohm's law is modified to include the Hall current so that 



ft> g T e 

~b7 



(JxB) 



E 



VxB + — V/? £ 
en P 



(2.7) 



in which co e is the cyclotron frequency, % e the electron collision time, cr the electrical con- 
ductivity, e the electron charge and p e the electron pressure. The ion-slip and thermoelec- 
tric effects are not included in (2.7). Further, it is assumed that co e T e ~ 01 and ©,T; <c 1, 
where <a, and T; are the cyclotron frequency and collision time for ions respectively. 

Inserting (2.1) in (2.5) and making use of (2.2), (2.3), (2.6) and (2.7) we obtain the 
following vector equation 



grad 



p|V| 2 + /?-ai ( V-V Z V + ^|A! 12 



-p[V f -Vx(Vx V)] 



: juV 2 V + «i [V 2 V r + V 2 (V x V) x V] 
+ (a, + a 2 ) div A 2 - ,' r ' - ^V. 



1 2 <r B§(l+iy) _ M, 
l + <p 2 K 



(2.8) 



in which V 2 is the Laplacian operator, V, = dV/dt, and |A] | is the usual norm of matrix 
Ai. If this model is required to be compatible with thermodynamics, then the material 
constants must meet the restrictions 17 81 



JU > 0, G5l > 0, (X\ + 0C2 = 0. 



(2.9) 



On the other hand, experimental results of tested fluids of second-grade showed that 
(X\ < and a\ + OC2 7^ which contradicts the above conditions and imply that such fluids 
are unstable. This controversy is discussed in detail in 1191 . However, in this paper we 
will discuss both cases, <X\ > and a,\ < 0. 

The velocity field for the problem under consideration is of the following form 



V(x,y,t) = [u(x,y,t), v(x,y,t), 0] 



(2.10) 



where u and v are velocity components in the x and y directions, respectively. 

Inserting (2.10) in (2.4) and (2.8) and making use of the assumption (2.9) we obtain the 
following equations 



du 
dx 

dp 

dx 



dv 

<5 = ' 



du 



(2.11) 



dt 



v(0 



li + at\ V 2 m - a.\vV 2 (o 

aBk(l+iq>) n 

-u u, 



l + (p 2 



K 



(2.12) 
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dp 



dt 



uco 



= ( pL + a\— J V 2 v + aiMV 2 o) 

ffB§(l + iq>) n 
l + <p 2 v ~k V ' 



(2.13) 



where (p — oo e T e is the Hall parameter and 



dv du 
dx dy ' 



ox 



mV 2 m + vW 



^ 'du dv\ 2 
5y <9x , 



(2.14) 
(2.15) 

(2.16) 



Remark 1. On setting a\ = 0,K — > °° and neglecting Hall effects in (2.12) and (2.13) 
we recover the equations for viscous fluid, on taking K — > °o and neglecting Hall effects 
we obtain the case [26 ] and on taking (X\ = and <p = in Hall effects we obtain the 
Brinkman model for porous medium. 

Equations (2.11 )-(2. 1 3) are three partial differential equations for three unknown func- 
tions u , v and p of the variables (x,y). Once the velocity field is determined, the pressure 
field (2.15) can be calculated by integrating (2.12) and (2.13). Note that the equation for 
the component w is identically zero. 

Eliminating pressure in (2.12) and (2.13), by applying the integrability condition 
d 2 p/dxdy = d 2 p/dydx, we get the compatibility equation 



da 



d d , 

U—+V— ) CO 

ox oy 



= ( ji + ai-^ ) V 2 co 



A 



ctB 2 (1 +iq>) 
l + <p 2 



dx 



0) 



-v— ) V-w 



(2.17) 



Let us consider the Stokes stream function: 



u = — , v- 
dy 



dy 

dx ' 



(2.18) 



where y(x,y) is the stream function. We see that the continuity equation (2.1 1) is satisfied 
identically and (2.18) in (2.17) yields the following equation 



-vV-{v,vV} 



>x + ai dt J v4 v / ~ ai [{^' v4 v / }] 

l + (p 2 r K r ' 



(2.19) 



in which 



V 2 -V 2 , © = -V 2 v/ 



(2.20) 
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and 



7 , , d Y d(V 2 Y) d W d(V 2 W ) 



ox oy oy ox 

is the Poisson bracket. 

Remark 2. The equation (2.19), for K — > °° and in the absence of Hall effects reduces 
to ESI. 

3. Solutions of some special types 

3.1 Solution of the type \j/(x,y) = %(x) + rj(.y) 

We consider the plane steady flow and examine the solution of (2.19) of the form 

V(x,y) = S(x) + Ti(y), (3.1.1) 

where % and 77 are arbitrary functions of the variables x and y respectively. Substituting 
(3.1.1) in (2.19) we obtain the following equation 

p fo'oorw - s W60] = mk w w + i lv wi 

+ a 1 [Tj'(y)| v (x)-rWr? V W] 
'aBg(l + i9) 



1 + <p 2 



(3.1.2) 



in which IV and V in the superscript indicates the fourth and fifth derivatives. 

We see that (3.1.2) is highly non-linear and its solution in the present form is not easy 
to obtain. In order to find its solution we assume the following 

%(x)=Ax+Be ax , (3.1.3) 

ri(y)=Cy+DeP (3.1.4) 

and obtain the following equation 

p + a 3 BCe ax - ab 3 BDe ax+by + a 3 bBDe" x+by ] 

-Ab 5 De b > + a 5 BCe ax 



■ H[a 4 Be ax + b 4 De by ] + ai 

'ff B§(l + iy) n 
l + (p 2 K 



-ab 5 BDe, ax+b ? + a 5 bBDe ax+b ? 
[a 2 Be ax + b 2 De by ], (3.1.5) 



where A,B,C,D,a and b are arbitrary constants. 

The following three equations are obtained from (3.1.5) 

paC = y.a 2 + a\a i C -H- ~, (3.1.6) 

K 

-pbA = ^b 2 -aib 3 A-H--, (3.1.7) 

K 

{b 2 -a 2 )[p-a l (a 2 + b 2 )]=0. (3.1.8) 
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From (3.1.6) and (3.1.7) we easily obtain the values of A and C, i.e. 
1 \ . H jx 



A = - 



p - (X\b 2 
1 



jib- 



b Kb 



p - a\a l 



H a 

[la 

a Ka 



(3.1.9) 
(3.1.10) 



where 



H = 



gBg(l+/«p) 

l + (p 2 

and (3.1.8) is satisfied if either 



b 2 -a 2 = 



(3.1.11) 



or 



p = ai(a 2 + b 2 ). (3.1.12) 

We have three different cases which we discuss separately as follows: 

Case I. b = a,p ^ a\(a 2 + b 2 ) 

The stream function given by (3.1.1), after using (3.1.3), (3.1.4), (3.1.9), and (3.1.10) 
becomes 



yr(x,y) 



(y-*) 

p - a\a 2 



\la 



a Ka 



-Be ax + De ay 



and from (2. 18) the velocity components take the following form 
1 



p - a\a L 



1 



p - a\a L 



H 




a 


Ka 


H 


M 


a 


Ka 



+ Dae ay , 
-Bae"*. 



(3.1.13) 

(3.1.14) 
(3.1.15) 



In order to find the pressure field (2.15) we substitute the velocity components (3.1.14) 
and (3.1.15) in (2.12) and (2.13) and then integrate the resulting equations to obtain 

p = PQ -pa 2 - LlBa 3 ye ax + (p - a x a 2 ) [a 2 Bye ax + a 2 DBe a{x+y ^] 

+ ai [B 2 a 4 e 2ax + DW ay - DBa A z a{x+y \ (3.1.16) 



where 



1 



jUfl 



H pL 



p — (X\a z |_' a Ka 
whereas the streamline for \ff = Cli is given by the following functional form 
-Be 



-xe + Q.\ 1 „ , 

y = Product Log 

e a 



a(--Be g *+jce+fli)/e 



(3.1.17) 
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where 



1 

1-Aa 2 



va- 



X 



v - 
Ka. 



jl/p as the kinematic viscosity, A = a,\/p is the second-grade parameter % = 
n N = cBq/p is the MHD parameter and <p = co e r e is the Hall param- 



with v 

N(l+i(p)/(l + (p 
eter. 

Streamlines are shown in figure 1 for B = D 
N = 0, co e r e = 0.1 and y/ = 15,20,25,30,40. 



-a = l, ju/p =0.5, ai/p = 0.1, # = 0.1, 



Case 2. = — a, p ^ 0£i (a 2 + fr 2 ) 
The expressions for y/,M, v, and /? are 



yr(x,y) = 



1 



p - a\a L 



H n 



jUfl — 

a Ka 



1 



M = 



p - aicr 
1 



p - aifl 2 



Ha 

a & 

jUfl 

a /fa 



(y + x)+Be ax + De- ay , 
-Dae- ay , 
Bat *, 



p = p Q -pa 2 - liBa 3 yt ax + (p - a x a 2 ) [a 2 Bye ax + a 2 DBe a ^} 
+ ai [B 2 a 4 e 2ax + D 2 a 4 e- 2ay - DBa 4 e a ^} 

and the functional form of streamline for \j/ — Q.2 is given by 

-Be ax -xe + Q. 2 , 1 



+ - Product Log 

a 



Da c a(Be ax +x£-n ? )/e 



(3.1.18) 
(3.1.19) 
(3.1.20) 

(3.1.21) 
(3.1.22) 



Streamlines are drawn in figure 2 for B = D = a = 1, fi/p = 0.5, a\/p = 0.1, K = 0.1, 
N = 0.5, co e T e = 1 and y/ = 15,20,25,30,40. 
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Figure 2. Streamline flow pattern for yf(x,y) = p _ l Uia i [m° _ f _ ks\ (>" + x ) 
+Be ax + De- a y. 



Case 3. b 2 -a 2 ^0 
We must have 

p = <X\(a 2 + b 2 ) 

and the expressions for \jf,u,v, and p are of the following form 



+ 



p - a\b 2 [ 

y 



p - a\a z 



H a 

Ha 

a Ka 



+ Be ax + De by , 



1 



p - a\a L 
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p - tt\b 2 



H a 

lla — 

a Ka 

u H ^ ' 
lib — 

^ b Kb 



+ Dbe by , 
-Bae a \ 



P = P0-^p[a 2 + b 2 }- fxBa 3 ye 



+ (p - a { a 2 ) [a 2 Bye ax + a 2 DBe ax+b y] 

+ cti [B 2 a A t 2ax + D Ve 26 ? - DBa A b 2 t ax+by ] , 



where 



(3.1.23) 
(3.1.24) 
(3.1.25) 



(3.1.26) 
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p - (X\a l 



\la 



H__y_ 

a Ka 



p - (X\b 2 



fib 



H jj_ 

~b~Kb 



whereas the functional form in this case for y/ = O3 is 

Be ax -x8-Q. 3 1 „ , T 

y = Product Log 

e b 



Db -HBe^-xS-a^/e 



, (3.1.27) 



where 



5 = 



1 



1-Ab 2 



vb-*-^- 
b Kb 
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Streamlines foiB=D = a = 1, ju/p =0.5,/? = -0.5, K = 2, 2.1, a x /p =0.5, -0J,N = 
0.5, OJgTe = 1 and y/ = 15,20,25,30,40 are depicted in figures 3, 4, whereas figure 5 is 
given for (X\/p = —0.5. 




1 2 3 4 5 

X 



Figure 3. y{x, y) = -xa + yb + Be ax + De hy . 
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The alternate forms of (3.1.23) may be written as 

H ju 



+ 



p - a\b 2 

y 



lib 



Kb 



a\b 2 

x 



H n 

[la 

a Ka 



-Be ax + De by , 



a\a L 



jib 



+ ■ 



y 



p - a,\a l 



H n ' 

~b~Kb 

H jj_ 

a Ka 



+ Be ax + De 



by 



(3.1.28) 



(3.1.29) 



Remark 3. The solution (3.1.16) with H = 0, K — > °° and tt\ = gives the Berker's solu- 
tion |4| and the Siddiqui's solutions 1261 can readily be recovered as a special case by 
taking H = and K — > °°. 

We now consider Riabounchinsky type flows in order to solve (2.15). 



3.2 Solution of the type \ff(x,y) =y^(x) 

In order to obtain another class of solution of (2.19) we substitute 

V(x,y)=yS(x) (3.2.1) 
into (2.19) and get the following equation 

p w" = ^ + ai u v - 4'4 iv ] -(#+£) r, (3.2.2) 

where 4 (x) is an arbitrary function of x, primes denote the derivative with respect to x. 
Integrating (3.2.2) once and equating the constant of integration equal to zero we obtain 

M r + p tor 2 - 44")] + «i \m IV - + 4" 2 )i - (h + 1) r = o. 



For the solution of the above equation we write 

4(.*) = S(l+Ae CT - v ) 



(3.2.3) 



(3.2.4) 



in which 8,<J and A are arbitrary real constants. Making use of (3.2.4) into (3.2.3) we 
have 



8 = 



1 



P -<J\G Z 

and thus from (3.2.1) 
\jf(x,y) 



— (h+—) 

a V KJ 



1 



p - a\o 



y(l + le ax ). 



(3.2.5) 



(3.2.6) 
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Figure 6. Streamline flow pattern for \jf(x,y) = p _„ i<T 2 [m ct ~~ F { H + §)] y(l+ 
Ae<"). 



The velocity components (2.18) and the pressure field (2.15) become 
1 



p - aio^ 



P - Gil C 



r 2„2' 



(l+Ae CT *), 
aXe ax . 



(3.2.7) 
(3.2.8) 



P = Pi +H°a (l - ^) Ae« - \p [a 2 (l A 2 e 2 -)] 
[^ + ^(3 + ^)^-], 



+ a\ 

where /?i is the reference pressure and 
1 



(3.2.9) 



p - (X\G l 



The streamline flow for y/ = Q.4 is given by the functional form 
i2 4 



(1 + Ae<™)e' 



(3.2.10) 



where 



1 



1-A<7 2 



VC7 



Figure 6 shows the streamlines for a = A = 1, ji/p = 0.5, <X\/p = 0.1, = 15, = 
0,y/ = 15,20,25,30,40. 

3.3 Solutions of the type y{x,y) — yt, (x) + 77 (x) 
Inserting 

\l/(x,y)=y${x) + r](x) (3.3.1) 
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in (2.19) we obtain the following equation 

-p\y(SV-K'") + &'!;"-&")] 

= ^ + 77 IV ) - a, 44 v ) + (77'4 IV - ^ v )] 

-(H + ^W + ti"). 

From the above equation, we have 

p[4'4"-4n + ^ IV -«i[r<? IV -^ v ]- (// + |)<?" = o 



(3.3.2) 



(3.3.3) 



and 



p fe'4" - 4t?"'] +MT? W - Oi [77'4 IV - 4»7 V ] - (# + I) l" = 0, (3.3.4) 



where £(x) and Tj(x) are arbitrary functions of its arguments. Integrating (3.3.3) and 
(3.3.4) and then taking the constants of integration equal to zero we have 

A i4'" + P[r 2 -^"]-ai[(-^ IV +24'4" , -4" 2 )]- 4' = o, 

(3.3.5) 

(ff + |)tj' = 0. (3.3.6) 



M Tj'" + p[7 ] '4'-47 ] "]-a 1 



4'T7'"- 477 IV 1 

+Tj'4'"-7]"4" 



We note that (3.3.5) is similar to (3.2.3). Its solution is given in (3.2.4). Substituting 
(3.2.4) into (3.3.6) we have 



aiA{\ + XQ ax )r\ + (fi - aiAXac ax )t]'" +A[{aia 2 - p)Xz ax - p)t)" 



(p - a x o 2 )AXoz ax - + 



(3.3.7) 



where 



A = 



1 



p - a x a z 



We note that it is not easy to obtain the general solution of (3.3.7). In order to find its 
solution we consider the following special cases: 

Case 1. OCi^0,(7 = l,A=0 
Equation (3.3.7) reduces to 



a x A x 77 IV + mt?'" - P A in " - (H + 1) t,' = 0. 



(3.3.8) 



We see that (3.3.8) is of fifth order and in order to solve it we reduce its order by putting 
rj' —A(x) such that (3.3.8) becomes 



aiAiA 1 " + juA" - pA x A' - (h + ^) A = 0. 



(3.3.9) 
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On substituting A(x) = P(x)e x , (3.3.9) takes the form 

Oi (3P' + 3P" + P" V + T" (2P' + P"Y -pP'e x = 0. 
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(3.3.10) 



Finally, P'{x) = R(x) converts (3.3.10) into a second-order differential equation 

CdR" + {}i/Ai+3a l )R' + {3a l -p+2[i/Ai)R = 0. (3.3.11) 
The solution of the above equation is 



R(x) = A 3 exp 



-c — Vc 2 -Ad 



x + A 4 exp j 



-c + Vc 2 -4d 



x, (3.3.12) 



where A3 and A4 are arbitrary constants and 

3aiAi+jU , (3<Xi- p)Ai + 2jx 
c— : , d = — 



aiAi 
1 

P-«l 



K 



«iAi 
H . 



In order to find 77 (x) we make backward substitutions and finally obtain the form 



7]{x) = 



A3 



+ ■ 



A 4 



m\{\+mi) m2(l+m2) 
where A, (i = 5,6) are constants of integration and 



e( 1+m2 ^+A 5 e x +A 6 x, (3.3.13) 



-Ad 



ni2 



-Ad 



From (3.2.4), (3.3.1) and (3.3.13) we get 



V(x,y) = 



y 



p-a\ 
+ 



M( 1- 

A 3 



1 



H 



K 

(\+m\)x ^_ 



+A 5 e x +A 6 
A 4 



(l+m 2 )x 



OTl(l+OTi) ni2(l+m2) 
The velocity components and pressure field are 



1 



p-«i 



Ml 



^l e ( 1+m i^ + ^e( 1+m 2> + A 5 e- r 
mi M2 



(3.3.14) 

(3.3.15) 
(3.3.16) 



P=P2-^p[A 1 ] 



+ 0£i 



3 e 2(l+mi> _|_ 2A 3 A 4 c (2+m,+nn)x _|_ A 4 e 2(l+m 2 )* _|_^2g2:r 



m\in2 



2A 3 A 5 (2+3m 1 +mf) (2+mi U i 2A 4 A 5 (2+3m 2 +m|) (2+»nk 
" r (2+mOmi c " r (2+m 2 )m 2 



(3.3.17) 



where /?2 is the reference pressure. 
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Figure 7. Streamline flow pattern for yr(x,y) = [jti (l — j) — -ff] +^56"* 

-1-4, 4. A 3 p (l+mi)s , A 4 (l+m 2 > 
+ A 6+ mi (l +mi ) e + m 2 (l+m 2 ) e 



The streamline for y/ = O5 is given by the functional form 

^3 „(l+m\)x 



where 



1 

£1 



+ mi {l+mrf e 

+ ,tl v> e' 1+m2 > + A 5 e* + A 6 x 

, m2(l+m2) 2 J ° . 



(3.3.18) 



1-A 



vM-- 



Streamline pattern is plotted in figure 7fora = A = l,/i/p = 0.5, a,\/p = 0.1, K = 0.5, 
N = 0, (p = 0.05,A3=A4=A5=A6=l, 1// = 15,20,25,30,40. 

Case 2. «i 7^ 0, a = 1, X ^ 
Equation (3.3.7) reduces to 

O1A1 (1 + Ae x )77 IV + (/i - aiAjAe^)^'" + A! [(01 - p)Xe x - p]rj" 



(p-a 1 )A l Xe K -(H+^ 



(3.3.19) 



To obtain the solution of (3.3.19) we try to reduce its order. For this purpose we put 
77' = A(x) which leaves (3.3.19) into a form which is one order less, that is 



01 (1 + Xe x )A"' + (ju/Ai - aiXe x )A" + [(ai - p)Xe x - p]A' 

A = 0. 



(p-aOAe*- + £ 



(3.3.20) 



Now substituting A(x) = P(x)e* in (3.3.20) and then P (x) = R(x) into the resulting 
expression, we get 

Kn(p-ai) 



ai(l+Xe x )R' 



K(n-H)-n 
2K/j,(p - a\) 



+ ai{3+2Xe x ) 



R' 



{2a x - p){\ + Xe x 



R. 



K{n-H)-n 

where we have taken the constant of integration equal to zero. 



(3.3.21) 



Inverse solutions for a second-grade fluid 
Subcase 1. The solution of (3.3.21) for X = is given by 

_ -Xi-^Xl-AX2 

R(x) = C\ exp 
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C2 exp 



-Xi + yJxf-4X 2 S 



(3.3.22) 



where C\ and C2 are arbitrary constants and 
Kjj,(p-ai) 



Xi = 



2Kii(p-ai) 3d- p 

3, -*2 = „,,.. ~ — + ■ 



ai[K(n-H)-n] ' ai[^T(ju-H)-ju] ai 
The backward substitution gives the value of rj (x) as 



(3.3.23) 



ri(x) = — 



Ci 



c 2 



mi(l+mi) m2(l+ffJ2) 
where C, (z = 1,2,3,4) are constants of integration and 



— e (1+m2) * + C 3 e* + C 4 , (3.3.24) 



Xi - Jxf-4X 2 -X { + Jxf-4X 2 
^ . m 2 = ~z • 



mi 



The stream function, the velocity components and the pressure field in this case are 
respectively given as 



y 



p-ai 

+ 



H 



m\(\ +mi) 



-C 3 e JC + C 4 

Ct_ 

m 2 (\ -\-Tn2) 



J}+m 2 )x 



p-ai 

- £± e (i+mi)* + — e (1+ ™ 2) * + C 3 e* 



mi 



P = /'3-2p^l] 



«2 



(3.3.25) 
(3.3.26) 
(3.3.27) 



+ai 



c ? c 2fl+mik 1 2ClC2 c (2+mi+m9);c | c % c 2(l+mi)x | ^2 C 

™2 



/n 1/712 



2v 



2CiC 3 (2+3mi+m 1 ) (2+mi)x , 2C 2 C 3 (2+3m 2 +m 2 ) (2+m 2 k 
1 (2+m! JthJ c ~ i (2+m 2 )m 2 



(3.3.28) 



where is the reference pressure and the streamline for = £2g is given by the functional 
form 



1 



Ci 



mi(l+mi) 2 



»(l+mi) 



_ e (l+m 2 > + C3e A : + C4 



(3.3.29) 



m 2 (l+m 2 ) 2 

Streamline pattern is plotted in figure 8 for a = X = l,ji/p = 0.5, 05i/p = 0. 1, K = 1.8, 
N= 1, (?>= l,Ci =C 2 =C 3 =C 4 = 1, ^= 15,20,25,30,40. 
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Figure 8. Streamline flow pattern for y/(x, y) = [jU ( 1 - ^ ) - H] + C?,e x + 
C 4 + -n\- , e( 1+ "'> + _ . e ( 1+ ^>. 

mi(l+mi) m2(l+f«2) 

Subcase 2. (General case). We now try to obtain the solution of (3.3.21) for Equa- 
tion (3.3.21) may be written as 

(1 + Xe x )R" + [X] + 2Xe x ]R' + [X 2 + X 3 Xe x ]R = 0. (3.3.30) 

In order to solve (3.3.30) we put 9 = e x to get the following equation 

(1 + X9)9 2 R" + (1 +Xi + 2X9)9R' + (X 2 +X 3 X9)R = 0, (3.3.31) 

where differentiation is with respect to 9. The solution for (3.3.31) is obtained through 
Mathematica and is given by 

R(9) = 9 X I 2 (- Xl -y/x 2 -4X 2 ^ 

C 5 9V X ?- 4X 22F1 { i*! , ±<& 2 , 1 + yjxf - 4X 2 , -OX } 

X 

+Q2F1 {1*3,3*4,1-^-4X2,-^} 

where 



*1 


= 2 


-Xi + y 


h 


-4X 2 -2^J\ 


-x 3 


<J>2 


= 2 


-Xi + y 


In 


-AX 2 + 2^f\ 


-x 3 


4> 3 


= 2 


-X\ — y 


M 


-4X 2 -2^/l 


-x, 



and 2^1 is the hypergeometric function defined in Appendix A. The stream function, 
velocity components and the pressure field can be obtained through the definitions of 

2 Fi[a,b,c,z\. 



Inverse solutions for a second-grade fluid 
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4. Concluding remarks 

In this paper, the exact solutions of non-linear equations governing the flow for a second- 
grade fluid in a porous medium are obtained by assuming different forms of the stream 
function (already used by various authors in different situations), in presence of a strong 
magnetic field. The expressions for velocity profile, streamline and pressure distribution 
are constructed in each case. Our result indicates that velocity, stream function and pres- 
sure are strongly dependent upon the material parameter a,\ of the second-grade fluid. It 
is shown through graphs that increase in (X\ leads to decrease in velocity and decrease in 
(X\ leads to increase in velocity (see figures 4 and 5). Also, the present analysis is more 
general and several results of various authors (as already mentioned in the text) can be 
recovered in the limiting cases. 



Appendix A 

Hypergeometric 2F1 [a,b,c,z] is the hypergeometric function 2F1 [a,b;c;z] and is the spe- 
cial case of the generalized hypergeometric function p F q [a;b;z] for p = 1 and q = 1. 
Hypergeometric function has the following properties: 

1. The 2F1 function has the series expansion 2F1 [a,b;c;z] — Y£=o ^(*)^* 

2. Hypergeometric 2Fi[a,b,c,z] has a branch cut discontinuity in the complex z-plane 
running from 1 to °°. 
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